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Abstract 

Background: Vector-borne disease transmission is dependent on the many nuances of the contact event between 
infectious and susceptible hosts. Virus acquisition from a viremic human to a susceptible mosquito is often 
assumed to be nearly perfect and almost always uniform across the infectious period. Dengue transmission models 
that have previously addressed variability in human to vector transmission dynamics do not account for the 
variation in infectiousness of a single individual, and subsequent infection of naive mosquitoes. Understanding the 
contribution of this variability in human infectiousness is especially important in the context of introduction events 
where an infected individual carries the virus into a population of competent vectors. Furthermore, it could affect 
the ability to detect an epidemic (and the timing of detection) following introduction. 

Methods: We constructed a stochastic, compartmental model to describe the heterogeneity of human viremia and 
calculate the probability of a successful introduction, taking into account the viremia level (and thus acquisition 
potential) of the index case on, and after, the day of introduction into a susceptible population and varying contact 
rates between the human and mosquito populations. We then compared the results of this model with those 
generated by a simpler model that has the same average infectiousness but only a single infectious class. 

Results: We found that the infectivity of the index case as well as the contact rate affected the probability of 
emergence, but that contact rate had the most significant effect. We also found that the interaction between 
contact rate and the infectiousness of the index case affected the time to detection relative to the peak of the 
epidemic curve. Additionally, when compared to our model that accounts for variable infectiousness, a model with 
a single infectious class underestimates the probability of emergence and transmission intensity. 

Conclusion: Understanding the interplay between individual human heterogeneity of infectiousness and the rate 
of contact with the vector population will be important when predicting the likelihood, detection, and magnitude 
of an outbreak. 

Keywords: Dengue, Transmission, Arbovirus, Heterogeneity of infectiousness, Human infectiousness, Dengue 
introduction, Stochastic modeling, Dengue model, Viremia 



Background 

Disease transmission is dependent on the many nuances 
of the contact event between infectious and susceptible 
hosts. In the dengue (DENV) transmission system, we can 
explicitly define the contact events as 1) transmission 
from Aedes mosquitoes with a productive, disseminated 
DENV infection (where virus is present in the saliva and 
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deposited upon mosquito probing and feeding); and 2) ac- 
quisition of virus by a naive Aedes mosquito from an in- 
fected individual with systemically circulating DENV [1]. 

Heterogeneity in the first of these- transmission from 
Aedes spp. to humans- can be further subdivided into the 
process of dissemination and, ultimately, successful trans- 
mission. Vector competence is the inherent ability of a 
mosquito to support replication of a virus for subsequent 
transmission and, as a static measure, is often calculated 
as the total proportion of DENV exposed mosquitoes that 
develop a disseminated infection capable of transmission 
at a particular time point. Vector competence has been 
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shown to vary based on viral and/or mosquito strains 
[2-7] as well as environmental and ecological conditions 
[2,4,8-11]. Further, there is variability of vector compe- 
tence at the level of the individual mosquito, and as such, 
vector competence is best evaluated as a dynamic process 
through a mosquito population [7,12], which takes into 
account the extrinsic incubation period (EIP, the time it 
takes for an arbovirus to replicate sufficiently through the 
vector for transmission to a vertebrate host) as well as the 
mortality rate of the mosquito [2,10,11,13,14]. 

The second of these contact events, the acquisition of 
virus from a viremic human to a susceptible mosquito is 
often assumed to be nearly perfect and almost always 
uniform across the infectious period. DENV transmission 
models have previously addressed variability in human to 
vector transmission dynamics mostly in the context of 
antibody-dependent enhancement (ADE). These studies 
have focused on secondary (rarely, tertiary or later) infec- 
tions and the theory that these infections have a higher in- 
tensity of transmission than primary infections due to 
ADE [15-19]. While these studies suggest that differential 
infectivity of DENV patients may be important, they are 
focused more on the average infectiousness of a group of 
individuals- severe versus not-severe DENV patients. 
These studies do not account for the variation in infec- 
tiousness of a single individual, and further, these studies 
do not directly address the consequences of this hetero- 
geneity: the infection of naive mosquitoes. Recently, a 
study in DENV infected patients showed that there is dif- 
ferential acquisition (and thus transmission) rates relative 
to individual host viremia among a cohort of symptomatic 
(both febrile and severe) individuals [20]. Understanding 
the contribution of this variability in human infectiousness 
is especially important in the context of introduction 
events where an infected individual carries the virus into a 
population of competent vectors. Since the introduction 
of DENV-1 into Florida in 2009-2010, the virus has be- 
come established in Key West and moved north into 
Miami and Martin county [21]. However, DENV has been 
detected in other locales in the United States such as Hous- 
ton and Brownsville, TX, as well as Suffolk County, NY 
[22-25] without long-term establishment of the virus. This 
is an indication that our understanding of the factors in- 
volved in successful emergence of DENV remains limited. 

While the factors affecting exposure to the virus are 
many (environment, seasonality, socioeconomic and life- 
style factors) [22,26], a critical component is the contact 
event between a viremic index case and the susceptible 
mosquito population. The key parameters affecting con- 
tact at the critical moment in introductory transmission 
into a completely susceptible population of mosquitoes 
are 1) the probability of acquisition of virus by a suscep- 
tible mosquito resulting from contact with the index case 
and 2) the contact rate between the mosquito and human 



populations. In addition, the interaction of these two 
transmission components defines the likelihood that an 
introduction will produce a fulminant, detectable outbreak. 

Herein, we investigated the role of dynamic human in- 
fectiousness in determining the emergence potential and 
subsequent detection of DENV-1 into a naive population 
of humans with the competent vector, Ae. aegypti. Add- 
itionally, while characterizing transmission explicitly in 
terms of the viremia curve and corresponding infectivity 
to mosquitoes, we determined how differential contact 
rates between the human and mosquito populations im- 
pacted the probability of detectable DENV-1 outbreaks. 

Methods 

We constructed a stochastic event-driven model of DENV- 
1 transmission that explicitly accounts for heterogeneity of 
human viremia, based on the compartmental framework 
of Susceptible-Exposed-Infectious-Recovered classes for 
humans and Susceptible-Exposed-Infectious classes for 
mosquitoes (Figure 1). Variation in human viremia is 
modeled by dividing the human infectious period into 
subclasses, each representing a single infectious day 
(ID) with an associated level of infectivity. Data from 
Nguyen et al. [20] (DENV-1) are used to parameterize 
infectivity levels on a per day basis. Specifically, we de- 
rive the probability of acquisition given a specific level 
of viremia and the course of viremia over several days 
of illness (DOI) by extrapolating curves to the discrete 
values obtained from Figure 2 and Additional file 1: 
Figure S4 in [20], respectively. These data were com- 
piled by allowing naive Ae. aegypti mosquitoes to feed 
on viremic individuals whose level of viremia was later 
determined via qRT-PCR. Mosquitoes were incubated 
and later tested for the presence of virus and thus suc- 
cessful transmission [20]. Based on these data, we di- 
vide the human infectious period into 10 days in the 
model (classes ID!-ID 10 ), with ID4-10 corresponding 
to the DOI2-8 from [20]. 

Characterization of human viremia and mosquito 
acquisition potential 

Our model assumes that viremia lasts roughly 7-10 days 
and viremia onset is, on average, 2 days prior to illness 
presentation, which is presumed to be DOI1 [20,27]. As 
the data from [20] spans DOI2-8, we wanted to extrapo- 
late the viremia values of DOI 1 and those two days prior 
to symptomology, with the assumption that three days 
prior to symptoms, the viremia was zero. We then define 
acquisition potential (q) as the probability that a naive 
mosquito feeds on a viremic individual and subsequently 
acquires a viral infection from that blood meal. To deter- 
mine the acquisition rates for the entirety of the infec- 
tious period, we did the following: 
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Figure 1 Model schematic. There were 10 human infectious classes each corresponding to a single infectious day (!D q _ qo ). Mosquitoes exposed 
to DENV-1 infectious individuals are infected at a rate q, associated with ID,. 



1. Fit a function to the viremia curve from [20] with 
the assumption of zero three days prior to symptom onset 
We determined that the shape of the viremia curve was 
similar to that of a Weibull probability distribution 
function, but obviously the scale was not. Through an 
iterative process of linear transformation of the data (to 
reconcile the scale of viremia vs. a Weibull function), 
we minimized the maximum difference between the 
transformed data and the Weibull curve. Details are 
available in the Additional file 2. Goodness of fit com- 
pared to the original data was verified using a boot- 



strapped version of the two-sample Kolmogorov-Smirnov 
test (ks.boot procedure, R package Match), where the null 
hypothesis states there is no difference in the probability 
densities of the two data vectors. The p-value indicated 
this was the case (p-value > .9999) and the fit was deter- 
mined to be sufficient. 

2. Estimated the logistic regression used to determine 
average acquisition given viremia level from Nguyen, et al 
Additional file 1: Figure S4 and extrapolate to ID 1-3. We 
determined that the logistic fit had parameter scale = 0.75 
and location = 6.5, and extrapolated values of acquisition 
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Figure 2 Extrapolated fit of probability of acquisition. The distribution of acquisition over infectious day (ID) given the relationship between 
infectious day and circulating viremia from [20]. Red dots are data points from [20], blue squares are the extrapolated points of acquisition for 
ID1-3 based on the curve fit (black line). 
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of 2%, 37% and 83% for ID1, ID2, and ID3 respectively 
(Figure 2), (fitting performed using R package Mass). An 
iterative process was performed to minimize the max- 
imum difference between the data and fitted values. 

DENV transmission model and experimental design 

To investigate the role of human infectiousness in the con- 
text of a potential introduction event, we initialized the 
model with a single index case. With respect to the viremia 
curve from [20], we varied the point at which the index 
case was introduced to the susceptible mosquito popula- 
tion. Thus we ran iterations of the model for each of the 
following scenarios. Introduction during: the latent period 
(IDO) and then in each of the infectious classes ID1-10. 

We simulated the resulting emergence (n = 2000 realiza- 
tions) for each parameterization of the model and calcu- 
lated the probability of successful emergence of DENV-1 
following the introduction via index case. We defined a 
successful introduction relative to the detection of cases, 
which we based on a recent sero-survey from Key West, 
FL when it was determined that likely only 2% of DENV-1 
cases presented clinically and were reported [21]. Thus, a 
successful outbreak may go undetected if the expected 
total number of DENV-1 cases was less than 50 (i.e. 2% < 
1). Because detection is a chance event, a 2% clinical pres- 
entation rate may actually result in more or less than 50 
cases. Here, to focus on the effects of heterogeneity of in- 
fectiousness and varying contact rates, we define a suc- 
cessful, detectable introduction resulting in emergence as 
one where there were at least 50 total cases resulting from 
the single index case. Given the recent introduction of 
DENV-1 in Martin County, FL in the summer of 2013 
[28], we based our total and susceptible populations on 
the two main urban centers of that county, Palm Springs 
and Stuart Island, which have a combined total population 
of approximately 38,000 people. The rest of the county is 
largely rural and thus less likely to be ecologically compat- 
ible with the DENV-1 vector Ae. aegypti 

Given that transmission from humans to mosquitoes is 
ultimately tied to contact between the two populations, we 
investigated the role of contact rate (a) on the probability of 
a detectable epidemic. Thus, we evaluated our model over 
a range of average contact rates: twice daily (a = 2), once 
daily (a = 1), once every 1.3 days (a = 0.75), once every two 
days (a = 0.5), approximately once every 3 days (a = 0.33), 
and once every 4 days (a = 0.25). 

Assumptions and parameterization 

Our model was run for a single year, without explicitly 
accounting for spatial or temporal variability in model 
parameters. Given the short time of simulation relative 
to human lifespan, birth and death rates of the human 
population were ignored and a constant population 
size was assumed. For female mosquitoes, we assumed a 



constant emergence rate and a constant per mosquito 
mortality rate, which leads to a stable female mosquito 
population. Since our efforts were performed with sub- 
tropical South Florida in mind, we assumed that the 
weather was permissive for Ae. aegypti year-long and 
that the abundance of Ae. aegypti did not differ signifi- 
cantly among seasons [29]. 

The model is defined by the events and corresponding 
transition rates in Table 1, and the parameter values (with 
sources) are given in Table 2. Stochastic realizations of the 
model were simulated using an algorithm that implements 
the tau-leap approximation to Gillespie s algorithm [30]. A 
time-step of 0.125 days was chosen to achieve computa- 
tional efficiency without sacrificing accuracy. All calcula- 
tions, curve fitting and model simulations were performed 
in R version 3.0.1. Stochastic simulation code is included 
as a Additional file 2 (SEI10R_stoch_model.R). 

Results 

Probability of detectable epidemic: relative infectiousness 
and contact rate 

We investigated how heterogeneity in infectivity during 
the course of a human infection affects detectable emer- 
gence potential of DENV-1 following introduction by an 
index case. Specifically, we varied the timing of introduc- 
tion relative to the infectiousness of the index case when 
s/he enters the susceptible population of humans and in 
the presence of a competent vector population. The 
probability of emergence was expectedly higher on those 
days when acquisition potential was highest, and was 
minimal for ID9-10. (It is important to note that these 
are not cumulative probabilities, but based on the cumula- 
tive effect of introduction on a particular ID. That is, for 
ID1, the probability reported is the probability of a suc- 
cessful introduction given that the index case comes into 
the susceptible population at ID1, but transmits at the ap- 
propriate rates from ID1 until cessation of viremia.) Prob- 
abilities of emergence are depicted in Figure 3 and given 
in Additional file 1: Table SI. 

A daily contact rate of 2 produced a detectable out- 
break 77% of the time if an index case was introduced 
into a naive population during the exposed period (ID = 
0). The probability remained relatively high (>70%) for 
ID 1-3, and then dropped to approximately 60%, 50%, 
36%, 21%, 7%, 2% and <1% for ID4-10, respectively. 

When the contact rate was once daily (a = 1), the risk 
of emergence was reduced to 49% for IDO and >40% for 
both ID 1-3. This probability did not drop off as steeply 
for ID4 and ID5, which had probabilities of 35% and 
27%, respectively. The probability of emergence was then 
19% for ID6, 10% for ID7, and 3.8% for ID8. Again, the 
probability of emergence was minimal (-1%) on ID9-10. 

The risk of emergence for a = 0.75 was also very similar 
for ID0-ID4 (39%, 37%, 36% 33%, and 28%, respectively) 
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Table 1 Definition of transition rates between compartments of the stochastic S-E-l-R model 



Event 


Change in state 


Transition rate 


Transmission from mosquito to human 


(S h/ E h )->(S h -1, E h + 1) 


aS h (l m /N h ) 


Onset of infectiousness in human 


(EhJDO-^CEh-1, IDt + 1) 


z e E h 


Transition from ID] to \D 2 


(ID 1; ID 2 )^(ID r 1, ID 2 +1) 


v ] \D ] 


Transition from ID 2 to ID 3 


(ID 2 ,ID 3 ) — > (ID 2 -1, D 3 + 1) 


v 2 ID 2 


Transition from ID 3 to ID 4 


(ID3JD4) — ► (ID 3 -1, ID4+I) 


V3ID3 


Transition from ID 4 to ID 5 


(ID4JD5)— ►(ID 4 -1, ID5+I) 


V4ID4 


Transition from ID 5 to ID 6 


(ID 5; ID 6 )^(ID 5 -1, ID 6 +1) 


V5ID5 


Transition from ID 6 to ID 7 


(ID 6 ,ID 7 )-»(ID 6 -1, ID7+I) 


v 6 ID 6 


Transition from ID 7 to ID 8 


(ID 7; ID 8 )^(ID 7 -1, ID 8 +1) 


v 7 ID 7 


Transition from ID 8 to ID 9 


(ID 8 ,ID 9 )-»(ID 8 -1, ID 9 +1) 


v 8 ID 8 


Transition from IDg to ID-iq 


(Do IDm) — > (IDo-1 IDm+11 


V9ID9 


Recovery in human 


(ID 10 ,Rh)^(ID 10 -1,Rh+1) 


V 10ID10 


Adult (female) mosquito recruitment 


(SJ->(S m + 1) 




Susceptible mosquito death 


(SJ-(S m -1) 


MS m 


Transmission from human to mosquito 


(S m ,E m )^(S m -1,E m + 1) 


aS m Ii(qilD/N h ) 


Exposed mosquito death 


(E m )-(E m -1) 


ME m 


Onset of infectiousness in mosquito 


(E m JJ->(E m -1,l m +1) 


bE m 


Infectious mosquito death 




ML 



Events and corresponding transition rates in the stochastic SEI 10 R-SEI model. For each event, we list only those states that change. Variable definitions are given in 
the Additional file 2 and parameter definitions and values are given in Table 2. 



Table 2 Values, definitions and sources for parameters used in modeling efforts 

Parameter (value) Definition 



Reference 



b" 1 (9 days) 

M" 1 (18 days) 

£ m (5000 females/week) 

z" 1 (4 days) 

(vmo)" 1 (1 day each) 

qi (.02)* 

q 2 (37)* 

q 3 (.83)* 

q 4 (-92) 

q 5 (-89) 

q 6 (-82) 

q 7 (-54) 

q 8 (-21) 

q 9 (-06) 

qio (-02) 



Contact rate 
Average extrinsic incubation period 

Average mosquito lifespan 
Emergence rate of adult mosquitoes 
Latent period (humans) 
Duration of each infectious subclass (ID1-10) 



Acquisition potential 



Varied (.25-2) 
[14] 

Approximated from [43-45] 
[31] 
[27] 
[20] 



*Extrapolated or directly from [20] 



Note that the emergence rate (e m ) results in an average mosquito density (per person 1 ) of <0.5. 

*Mosquito acquisition (q^) values were extrapolated by curve fitting while q 4 . 10 were taken directly from Nyguen, et al. [20] as described. 
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Figure 3 Probability of successful, detectable emergence of 
DENV-1. Probability of DENV-1 emergence (y-axis) given the 
introduction via an index case on infectious day (ID, x-axis) for several 
contact rates (a) (defined by the color bar on the right). Error bars 

correspond to 95% binomial confidence intervals: p ± 1 .96 ^/p(1 -p)rr 1 . 
\ ) 



and then dropped to 22%, 14%, 8%, and 3% for ID4-8, and 
then became <1% on ID9 and ID 10. The probability of 
emergence when contact rate was once every two days 
(a = 0.5) was (in order of ID0-ID10) 24%, 24.5%, 23%, 17%, 
12%, 8%, 4%, 1.6% and for ID9-10, <1%. For all days but 
IDO, the probability of emergence when a = 0.33 was con- 
sistently below 15% (range = [<1-12%]). For a = 0.25, the 
probability of emergence was consistently below 5% 
(range = [.1-4.5%]), for IDO-9 and 0 for ID10. For contact 
rates less than a = 0.5, transmission was either very tenu- 
ous (a = 0.25) or did not peak for all simulations within 
365 days (a = 0.33). 

Examination of the relationship between the probabil- 
ity of successful emergence and contact rate found that 
it followed an approximately negative exponential curve 
(Additional file 1: Figure SI). 

Window of Infectiousness 

We were interested in determining the extent to which 
the length of an index case s productive infectivity was 
altered due to differential contact rate. That is, at what 
contact rates are periods of especially low viremia able 
to result in a successful transmission event (from a hu- 
man to a mosquito) and, conversely, at what point does 
the contact rate become such that only peak days of 
viremia contribute to the detectable emergence potential 
of DENV. The rate of contact between the index case 
and susceptible mosquito population altered the "win- 
dow of infectiousness" of that index case. Indeed, the 
probability of emergence during the days of lowest ac- 
quisition potential (ID 8-10) was consistently less than 
10%. For ID9 and ID 10, the probability was even lower: 



at its highest only 2.4%, and otherwise <1.1%, meaning 
that periods of similarly low viremia do not contribute 
meaningfully to the emergence potential of DENV-1. In 
addition, with a contact rate of a = 0.25 the window of 
infectiousness was truncated by two days as ID9 and 
ID 10 were non-productive and did not result in a suc- 
cessful emergence event. 

Proportion of population infected 

The magnitude of resulting outbreaks during the time 
over which the simulations were realized (365 days) was 
also examined. When contact was equal to or greater 
than a = 0.5, there was no apparent effect of ID on the 
cumulative proportion of the population that became in- 
fected during those 365 days. For a = [0.75, 2], there was 
saturation of the population (>99%) (Additional file 1: 
Table SI). That is, almost all susceptible individuals be- 
came infected, assuming no other forces interrupted 
transmission (such as environmental phenomenon, etc.). 
At a contact rate of a = 0.5, the proportion of the popu- 
lation infected was approximately 97%, regardless of the 
day of introduction. At a = 0.33, the proportion of the 
population infected ranged from 47-67.5%; and at a = 
0.25, the range was 1.5-.3.7%. 

Time to peak cases and transmission intensity 

We were interested in the time between probable detec- 
tion of a DENV-1 outbreak and the epidemic peak. To 
examine this, the peak number of cases was centered at 
time = 0 for each model realization that met the criteria 
for an outbreak (cumulative cases >50), and then the 
average number of cases for each time point relative to 
time = 0 was calculated (time = time actua i- time cen tered)- 

We assumed that there needed to be 50 cases before a 
clinical case would be detected, given the subclinical rate 
we extrapolated from the emergence of DENV-1 in Key 
West [21]. The time to peak cases (relative to detection) 
is then defined as the time between peak number of 
cases (at time cen tered) and likely detection of the first 
clinical case. We determined that contact rate was 
understandably tied to the average time to peak cases, 
but also that the variance in time to peak cases due to 
ID was related to contact rate. When the contact rate 
was greater than a = 0.25, outbreaks were well defined 
and the time between detection and peak was relatively 
stable (Figure 4, Additional file 1: Table SI). However, 
when a = 0.25, epidemics were erratic and thus this 
metric of detection relative to peak timing was deemed 
uninformative (Figure 4, Additional file 1: Table SI). 

A shorter time between detection and peak indicates 
either an epidemic where the preliminary outbreak is 
more intense in nature (versus a slow burn-in), or else is 
an artifact of an already tenuous transmission system. 
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To demonstrate this, we also investigated the cumula- 
tive proportion of the population that was infected at (or 
before) the peak of the epidemic. At the highest contact 
rate, the time to peak relative to detection was 4 weeks 
and the proportion of the population infected was be- 
tween 65 and 69%. When the timing of the peak relative 
to detection was still relatively short (7-8 weeks for a = 1 
and a = 0.75, respectively), the average proportion in- 
fected was approximately 60% at peak (range 58.6- 
64.1%). When the time to peak was between 15 and 
16 weeks (a = 0.5), the proportion infected at peak de- 
creased to approximately 52%. And at a = 0.33, when 
peak timing was between 35 and 37 weeks, the propor- 
tion infected at peak was between 33 and 40%. The pro- 
portion infected at peak for a = 0.25 ranged between 1 
and 4%, further illustrating the tenuous nature of trans- 
mission at this contact rate. However, for the other con- 
tact rates, the inverse relationship between time to peak 
(relative to detection) and the average proportion of the 
population that had been infected suggests there is a 
relationship between time to peak from detection and 
transmission intensity. 

Comparison with simpler SEi 1 R model 

We compared the results of our SEI 10 Rmodel structure 
to those of a more traditional SE^R model, which has a 
single infectious class with an average duration of 
10 days. We parameterized acquisition as the average of 



the rates for each of the 10 days in our multi-class 
model (q =0.468). This was done to enable a more direct 
comparison of the models because the reproductive num- 
ber of each model would be identical (see Additional files 
1 and 2). Results of the simpler model are limited to pre- 
dictions based on the exposed and single infectious period 
(equivalent to ID0 and ID1 above). Because introductions 
on ID0 and ID1 gave similar results, we focus on com- 
parisons between the introduction events occurring 
during the latent period (ID0). When considering the 
probability of emergence, the simple model predic- 
tions were consistently below those of the SEI 10 R for 
contact rates higher than a = 0.25, and the difference 
in emergence probability ranged from 2.25% (a = 0.33) 
up to 16.3% (a = 2). For a = 0.25, the simple model overesti- 
mated the probability of emergence by 0.15%, but this is 
again likely an artifact of tenuous transmission over the 
time period of simulation (365 days). Results from this 
model are given in Additional file 1: Table S2, Figures S4 
and S5. 

Although the two models had near identical results 
when estimating the cumulative proportion of infected 
individuals, there were also differences in the timing of 
the peak relative to detection. The simple SE^R model 
predicted a longer time to epidemic peak (relative to de- 
tection) when contact rates were sufficient such that the 
epidemic had always peaked during the 365 days: the 
disparity in timing estimates differed by +1 week (a = 2,1), 
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+2 weeks (a = 075), and +3 weeks (a = 0.5). The SE^R 
model also predicted time to peak relative to detection 
would be longer by one week when a = 0.33, though the 
outbreaks at this contact rate did not all peak within 
365 days and such a temporal metric is likely skewed in 
this case. The differences in predictions between the two 
model formulations could have significant implications 
for control strategies, which are usually reactionary, but 
carefully considered due to resource constraints. 

Discussion 

Our results highlight the importance of accounting for 
heterogeneity at the interface of human to vector DENV 
transmission. Our model framework has the granularity 
of day-to-day contribution of variable human infectious- 
ness and quantifies how the infectious day of an intro- 
duced case can impact outbreak potential. Indeed, there 
have been many failed DENV introductions into South 
Florida [28], which may be due, in part, to travelers being 
in the latter stages of viremia. While virus was detected in 
these travelers [28], our model illustrates a possible reason 
as to why these individuals did not produce viable and 
persistent chains of transmission. In addition, accounting 
for heterogeneity in human infectiousness leads to more 
accurate predictions of emergence potential and peak tim- 
ing compared to a model assuming only a single infectious 
class with the same average acquisition rate. In particular, 
for high contact rates, the simple model consistently 
underestimated the emergence potential and overesti- 
mated the time to peak cases from detection. 

The magnitude of the effect of variation in human infec- 
tiousness is significantly modified by human-vector con- 
tact rates. Differences in contact rates can be indicative of 
several broad factors: mosquito behavior differences, eco- 
logical and environmental stressors, and mosquito avoid- 
ance behavior by people [32]. Given the distribution of the 
primary vector (A e. aegypti) as well as a competent alter- 
native vector (Ae. albopictus), not only in tropical and 
sub-tropical regions but in urban centers like NYC and 
Atlanta [33], reducing human-mosquito contact is likely 
the only feasible means of transmission interruption. This 
can be achieved by several means: 1) insecticide-based 
vector control; 2) breeding source reduction; and 3) im- 
proved infrastructure in areas of lower socio-economic 
status (e.g., window/door screens). The recent DENV-1 
emergence in Martin County was in a more established 
neighborhood where people were anecdotally reported to 
have increased contact with the mosquito population due 
to extended periods spent outdoors and the absence of in- 
frastructure (such as screens and air conditioning) [34]. 

In fact, Martin County, FL had 23 cases of autochthonous 
DENV-1 transmission, 15 of which were residents of the 
county, in the summer of 2013, beginning in July and last- 
ing through late September (approximately 10-12 weeks) 



[28]. These cases were all reported to the Florida Depart- 
ment of Health, and were clinically presenting individuals. 
Assuming that the rate of clinical presentation of DENV-1 
in Martin county would be similar to that of Key West 
following DENV-1 introduction there (-2%) [21], then 
Figure 5 shows the extrapolated upper limit of total (cumu- 
lative) cases (clinical + subclinical = 1100 cases ~ 2.9% of the 
likely at-risk population). Our simulations suggest that, 
given the relative environmental stability of the area for the 
duration of the epidemic and forgiving other large-scale ex- 
trinsic forces, the average daily contact rate across the at- 
risk population (Palm City and Stuart Island) is more likely 
to have been closer to 0.25 than higher values. However, as 
DENV is often associated with clusters, it is likely that 1) 
the spatial scale of the epidemic was more limited, 2) the 
at-risk population was smaller and 3) the contact rate likely 
higher within that population [34-36]. 

Caveats 

There are a few model assumptions that deserve further 
consideration. First, we consider only the primary DENV 
vector, Ae. Aegypti, although Ae. albopictus is a competent 
vector for DENV, feeds on humans, and is present in the 
United States [29,37]. A major motivation of this was that 
the data on which the acquisition fit and subsequent mod- 
eling was performed comes from a human-Ae. aegypti ex- 
periment [20]. While no direct acquisition studies involving 
Ae. albopictus have been published, there are important dif- 
ferences in the vector competence of this secondary vector 
[38]. However, parameterization of this model could easily 
accommodate an effort focused on Ae. albopictus with ap- 
propriate attribution given to foraging preference, contact 
rates, and ecological niche of this species [39]. In terms of a 
two-mosquito model, that is the focus of on-going work in 
our laboratory and out of the scope of this particular effort. 

Second, our assumption regarding presentation and 
detection rate (2%) is based on the Key West, FL experi- 
ence [21] because Key West represents the most likely 
scenario to be repeated on subsequent introduction: 1) a 
mostly naive population, so the risk of secondary DENV 
is less than endemic areas, 2) relatedly, the clinical pres- 
entation rate (associated with secondary infections) will 
be more similar than in DENV-endemic areas and 3) en- 
demic areas are often set up for active disease surveil- 
lance and the medical community is trained to look for 
DENV symptoms. As DENV is a newly emergent patho- 
gen in the United States, our surveillance and reporting 
infrastructure is still developing and thus best based on 
the most recent event: Key West, FL in 2009. 

Finally, future consideration should be given to viremia 
differences within and among serotypes of DENV, as evi- 
denced by the experimental data in Nguyen et al. [20]. 
Another study has also reported differences among sero- 
types as well as an approximate log difference in DENV-1 
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Figure 5 Recent DENV-1 emergence in South Florida. Map of Florida showing Key West where DENV-1 was introduced in 2009-2010 and Martin 
County, which experienced a DENV-1 introduction and emergence in summer 2013. The bar graph depicts the number of reported cases from Martin 
County and the line graph depicts the extrapolated total proportion of residents infected using the sero-conversation rate of Key West. 



average viremia between primary and secondary infections 
[40]. While secondary infections would not be relevant for 
emergence of DENV, it would be very relevant for predic- 
tion of the incursion of a second serotype into an area that 
has only experienced one. This is likely going to be rele- 
vant in the near future as DENV-1 has been involved in 
autochthonous transmission in at least two counties in 
Florida, but other serotypes are often introduced by trav- 
elers [28], 

Conclusions 

DENV has been purported to be a pathogen spread pri- 
marily to new areas by people, rather than its vector, 
which has a relatively short flight range [41-43]. The 



relative infectiousness of an individual is especially im- 
portant when that individual is a potential index case, 
introducing the pathogen to a new area. In the case of 
DENV and other vector-borne diseases, the important 
parameters are the infectiousness of the viremic individual 
to a naive, competent vector and the contact between 
the infectious individual and a competent, sufficiently 
large mosquito population. In this paper, to quantify the 
emergence potential of vector-borne pathogens such as 
DENV-1 in predominantly susceptible populations, we 
incorporate a critical parameter (acquisition potential) 
that accounts for the viremic state of the index case. 
Further, by assessing differential contact, we address the 
role of other factors (like the suitability of the ecology 
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that such an index case might inhabit or arrive in) in facili- 
tating transmission and emergence. We observe that the 
probability of a successful introduction is greatly affected 
by the timing of introduction relative to the infectiousness 
of the index case and that when this heterogeneity is ig- 
nored, estimates of emergence probabilities are under- 
estimated. Further, the contact rate between the human 
and mosquito populations contributes to the emergence 
potential and the intensity of transmission. In general, 
our results highlight the importance of data-driven 
parameterization of acquisition based on the dynamic 
viremia process, and accounting for transmission het- 
erogeneity in terms of contact. 
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